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H I 21 cm opacity fluctuations power spectra towards Cassiopeia A 



o 

(N 

(N 
(N 

< 

o 

43 

Oh 
6 

(N 
> 
(N 
Os 

(N 
O 
O 



Nirupam Roy u , Jayaram N. Chengalur 2 *, Prasun Dutta 3 * and Somnath Bharadwaj 3i 

1 National Radio Astronomy Observatory, 1003 Lopezville Road, Socorro, NM 87801, USA 
^National Centre for Radio Astrophysics, TIFR, Post Bag 3, Ganeshkhind, Pune4Il 007, India 

3 Department of Physics and Meteorology & Centre for Theoretical Studies, IIT Kharagpur, Kharagpur 721 302, India 



Accepted yyyy month dd. Received yyyy month dd; in original form yyyy month dd 



ABSTRACT 

The angular power spectrum of H 1 2 1 cm opacity fluctuations is a useful statistic for quantify- 
ing the observed opacity fluctuations as well as for comparing these with theoretical models. 
We present here the H I 21 cm opacity fluctuation power spectrum towards the supernova 
remnant Cas A from interferometric data with spacial resolution of 5" and spectral resolution 
of 0.4 km s _1 . The power spectrum has been estimated using a simple but robust visibility 
based technique. We find that the power spectrum is well fit by a power law P T {U) = U a 
with a power law index of a ~ —2.86 ± 0.10 (3a error) over the scales of 0.07 — 2.3 pc for 
the gas in the Perseus spiral arm and 0.002 — 0.07 pc (480 — 15730 au) for that in the Local 
arm. This estimated power law index is consistent with earlier observational results based on 
both H I emission over larger scales and absorption studies over a similar range of scales. 
We do not detect any statistically significant change in the power law index with the velocity 
width of the frequency channels. This constrains the power law index of the velocity structure 
function to be f3 = 0.2 ± 0.6 (3er error). 
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1 INTRODUCTION 

It is now well established that the H I 21 cm opacity in the Galaxy 
shows measurable small scale structure. Very Long Baseline Inter- 
ferometry (VLBI) observations show opacity fluctuations on an- 
gular scales as small as 20 milli arc seconds corresponding to a 
projected separation ~ 10 au , if one assumes tha t the opacity vari- 
ations occur in a thin screen terogan et al.| [2005). Similarly, multi- 
epoch observations of some high velocity pulsars sho w evidence 
for H I 21 cm opacity variation on scales of 5 - 100 au jFrail et all 
1 1994) . Several other pulsa rs however do not show opacity varia- 
tions Jjohnston etail f2003). leading to suggestions that the occur- 
rence of fine scale structure in the Galaxy may be rare. Much of 
these observed fluctuations on au scales were earlier believed to 
originate in H I "clouds" with densities ~ 10 4 — 10 s cm -3 . It 
is hard to explain the existence of such structures in equilibrium 
with other, orders of magnitude lo wer density, compo nents of the 
diffuse interstellar medium (ISM). iDeshpandd (2000) has shown 
that the observed small scale transverse variations have contribu- 
tions from structure on all scales and that the fluctuations at larger 
spatial scales could, in projection, produce the measured small 
scale fluctuations. Conversely, there are other observational and 
numerical simulation results supporting the existence of tiny scale 
structures which have a size of ~ 3000 au, with a density of ~ 
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100 cm" 3 (e.g. [ Braun & Ka nekar 2005; Vazque z-Semadeni et al.l 
2006; Hennebelle & Audit 2007). With an evaporation timescale of 
~ 1 Myr, these structures can survive only if the ambient pressure 
is much higher than the normal value or if they are formed continu- 
ously on a comparable timescale. Details regarding their origin and 
physical prop erties are, however, still not understood. On somewhat 
larger scales. iDeshpande et al.l ( 120001) measured opacity variations 
of the H I absorption across Cassiopeia A and Cygnus A, and found 
that the power spectrum of these fluctuations could be fairly well 
represented by a power law. The slope of the power law was sim- 
ilar for the absorption in the Perseus arm (towards Cas A) and the 
Outer arm (towards Cygnus A), but substantially different for that 
from absorption towards Cygnus A arising in the Local arm. 

The scale free behavior of fluctuation power spectra of a 
variety of tracers (H I 21 cm emission intensity, dust emission) 
is the main observational evidence for the existence of turbu- 
lence in the atomic ISM. The slope of the power law of H I 
21 cm emission and absorption in our own Galaxy, H I 21 cm 
emission from the Large Magellanic Cloud , the Small Magel- 
lanic Cloud and DDO 210 are all ~ —3 iCrovisier & Dickey 
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1983l:lGreeril993l:IStanimirovic et al.l i999];lElmegreen et al.l200ll : 
Begum et alj|2006t) . Recently iDutta et alj J2009allbl) have reported 



that the H I 21 cm intensity fluctuation power spectra for a sam- 
ple of dwarf and nearby spiral galaxies have a power law index 
of ~ —2.6 and ~ —1.5 on scales respectively smaller and larger 
than the scale height of the galaxy disk. This is interpreted as the 
effect of a transition from three dimensional to two dimensional 
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turbulence at scales larger than the disk thickness. The observed 
fractal structure of H I in several dwarf galaxies in the M 8 1 group 
is also consistent with the self-similar hierarchical structure of the 
turbu lent ISM without any preferred length scale (Westpf ahl et al.l 
1 19991) . 

On the theoretical side, fine scale structure in the neutral gas is 
naturally expected in turbulent models of the ISM. The fluctuations 
in the H I 21 cm opacity in a particular velocity range depend on 
the fluctuations in the densit y, spin temperature and velocity of the 
gas. lDeshpande et"al] fcood) show that in the case of small fluctu- 
ations of density, the dependence on temperature is small. Turbu- 
lence however, gives rise to fluctuations in both the density and ve- 
locity of the gas and both of these cont r ibute to the observed opac- 
ity fluctuations. lLazarian & Pogosvanl j200(l) show that the slope 
of the observed power spectrum changes depending on whether the 
H I 21 cm emission is averaged over a velocity range that is large 
("thick slices") or small ("thin slices") compared to the turbulent 
velocity dispersion. Observations with high enough spectral reso- 
lution can hence disentangle the statistical properties of the velocity 
and the density fluctuations. 

Cassiopeia A or Cas A (Gil 1.7—2.1) is in the constellation 
Cassiopeia in the Galactic second quadrant. It is a shell type su- 
perno va remnant of di ameter 5' at a distance of approximately 3.4 
kpc dReed et alll 19951) . and is one of the brightest radio sources 
in the sky (flux density of ~ 2720 Jy at 1 GHz). The radio im- 
ages show it to have a clear shell like structure along with com- 
pact emission knots, the shell thickness estimated from the ra- 
dial brightness profile is ~ 30". There have been many stud- 
ies of the ISM towards Cas A using H I 21 cm absorption (e.g. 
[ Clark et alJI 19621; Iciarkl 1 1 9651 : ISchwarz et ai]|l986r . iBieging et al] 
1991; Revnoso et al. 1997) and a variety of other trace rs like H2CO 
and OH dGoss et al]|l984( : lBieging & Crutchen 1 ! 19861) all of which 
reveal small scale structure. A recently developed formalism for 
statistically robust extraction of the fluctuation power spectrum 
from interferometric data is used here to re-examine the issue and 
to constrain the H I 21 cm opacity fluctuation power spectrum to- 
wards Cas A. This method is particularly useful because one can 
completely avoid the complications of deconvolution and imaging 
the strong background source. We note that the op acity power spec- 
trum t owards Cas A has been obtained earlier by Des hpande et al.l 
(2000). They used the data from the Very Large Array (VLA) B, 
C and D configurations but with only 18 antennas in each configu- 
ration due to software restrictions. Estimating the power spectrum 
using the visibility based formalism for the same line of sight will 
allow us to compare the results from two very different methods to 
cross-check different techniques. Section $2] briefly describes the 
data while our methods of estimating the intensity and opacity fluc- 
tuation power spectra are outlined in Section Sj3]and ^respectively. 
Section ^contains the results and discussion, and we present con- 
clusions in Section 36] 



2 DATA 

The Giant Metrewave Radio Telescope (GMRT) L-band (21 cm) 
receiver was used to observe the H I in absorption towards super- 
nova remnant Cas A. The observation was carried out on Decem- 
ber 03, 2006 and the total duration was about 13 hours with on- 
source time of about 8 hours. VLA calibrator sources 2148+611 
and 2355+498 were used for phase calibration. One of the phase 
calibrators was observed for 3 minutes for every 20 minutes obser- 
vation of Cas A. Standard flux calibrators 3C48 and 3C286 were 
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Figure 1. Integrated H I 2 1 cm absorption spectrum towards the supernova 
remnant Cas A. The continuum corresponds to the total flux density of 
Cas A. The wide feature at negative LSR velocity is the absorption pro- 
duced by gas from the Perseus spiral arm and the absorption near zero LSR 
velocity is from the Local arm gas. Power spectra for channels marked as 
A, B, C and D are shown in Figure|2] 



observed for about 20 minutes in every 4 hours during the observa- 
tion. Frequency switching was used to calibrate the spectral base- 
line — the flux calibrator scans (which were also used for bandpass 
calibration) were observed at frequency offset by 5 MHz from that 
of the absorption line. A total baseband bandwidth of 0.5 MHz di- 
vided into 256 frequency channels centered at 1420.5973 MHz was 
used for the observation. This corresponds to a velocity resolution 
of ~ 0.4 km s _1 per channel and the total bandwidth was suffi- 
cient to cover the H I 21 cm absorption produced by gas both from 
the Perseus arm and the Local arm. Data analysis was carried out 
using standard AIPS. After flagging out bad data, the flux density 
scale, instrumental phase and frequency response were calibrated. 
The calibrated visibility data of Cas A were then used to estimate 
the angular power spectra and the errors. The integrated H I 21 cm 
absorption spectrum towards Cas A derived from the GMRT obser- 
vation is shown in Figure (T). The wide feature at negative velocity 
with respect to the Local Standard of Rest (LSR) is the absorption 
produced by gas from the Perseus spiral arm and the absorption 
near zero LSR velocity is from the Local arm gas. 



3 THE INTENSITY FLUCTUATION POWER SPECTRUM 

The angular power spectrum Ps{U) of the sky brightness fluctua- 
tions SS(l, m), is the Fourier transform of the autocorrelation func- 
tion £s 



P S (U) =Ps(u,v) = / Us{l,ra)e- 2 ^ ul+vm) dldm. 



(1) 



Here (I, ra) refer to the directions on the sky, (u, v) 
the inverse angular separation and £s is given by 

£ s (l - I', m - m') = (SS(l, m)SS(l',m')} . 



U refer to 



(2) 



The angular bracket refers to an ensemble average over different 
realisations of 8S(l,m) which is assumed to be stochastic in na- 
ture. The fluctuations are also assumed to be statistically isotropic 
which implies that Ps{U) depends only on U —\ U \. 

The visibility V V (U) measured in radio interferometric obser- 
vations directly probes a Fourier mode of the sky brightness flue- 
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tuation 8S(l, m). Here the inverse angular scale U directly corre- 
sponds to a baseline, the antenna separation projected in the plane 
perpendicular to the direction of observation, measured in units 
of the observing wavelength A. It is common practice to express 
the dimensionless quantity U in units of kilo wavelength (kA). The 
Fourier relation between V„ (Z7) and 5S(l, m) allows us to identify 
each baseline with an inverse angular scale, and to directly estimate 
the angular power spectrum from the measured visibilities. 

Detailed discussions of the visibility based power spec- 
tral e stimat or that we use he re can be found in Begu m et al.l 
fcOOd^ and TPutta et alj fe009ah . Briefly, the estimator P(f7) = 
{ V V (U)VZ(U + AU) ), is obtained by correlating every visibility 
V V (U) with all other nearby visibilities V„(U + AU), i.e. those 
which lie within the disk defined by \ AU\ < (ttOo) -1 , where do is 
the angular radius of the source. In the case of Cas A, (9o ~ 2.5', 
i.e. the angular radius of the source. Since the fluctuations are as- 
sumed to be statistically isotropic, the correlations are averaged 
over different U directions to obtain the final estimator. To increase 
the sig nal to nois e ratio w e further average the estimator in bins 
of U. iDutta et all {2009a) show that the expectation value of the 
estimator P(U) is the convolution of the power spectrum P(U) 
with a window function |W^,(J7)| . The window function, which 
quantifies the effect of the finite angular extent of the source, is 
peaked around (7 = and has a width of order (nOo)^ 1 beyond 
which IW^fV)! 2 ~ 0. For the optically thin shell-type geometry 
of Cas A the convolution does not effect th e shape of the power 
spectrum beyond a baseline U m ~ 1.6 kA l IRov et alj|2009l ), and 
we may directly interpret the real part of the estimator P(U) as 
the power spectrum P(U). The estimator P(U) also has a small 
imaginary component arising mainly from noise. We compute the 
la error bars for the estimated power spectrum by assuming the 
error to be the quadrature sum of contributions from two sources 
of uncertainty, viz. the sample variance and the system noise. At 
small U the uncertainty is dominated by the sample variance which 
comes from the fact that we have a finite and limited number of 
independent estimates of the true power spectrum. At large U, it is 
dominated by the system noise in each visibility. 

The longest baseline in our observation (~ 25 km) corre- 
sponds to an angular resolution of ~ 2.5", but our power spectrum 
estimation is restricted to U ^ 50 kA (angular scale ~ 5") be- 
cause the larger baselines do not have an adequate signal to noise 
ratio. Similarly for U ^ 1.6 kA (angular scale ~ 2.5') the sample 
variance leads to large uncertainties in the power spectral estimate. 
Note that since the estimator is local in U, the estimate of the power 
spectrum over this specific range in U is not affected by the miss- 
ing zero spacing baselines or the cut-off at large U. Certainly these 
cut-offs and consequent Gibbs phenomena will produce significant 
ringing and distortion in the image. However, since the power spec- 
trum is estimated directly in the uv plane, our method, unlike image 
based techniques, has completely avoided these problems. 

Figure[2]shows the estimated power spectrum for four differ- 
ent channels with H I absorption features of different optical depth. 
Note that since channel A has no H I absorption, it just shows inten- 
sity fluctuations power spectrum of the background source Cas A. 
On the other hand, channel D has inadequate signal to noise to 
constrain the power spectrum due to very large optical depth. The 
power spectrum Pj c of the continuum emission from Cas A, seen in 
the ch annels without H I absorption, has been studied in lRov et al.1 
(2009). This power spectrum is well fit by a broken power law 

Pi c = CU~ 222 for U < 10.6 kA 

= 10.6 101 CT/~ 3 ' 23 for U > 10.6 kA (3) 
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Figure 2. Power spectra for channels marked as A, B, C and D in Figure[TJ 
Typical error bars are plotted for only one channel for clarity. Since channel 
A has no H I absorption, it just shows intensity fluctuations power spectrum 
of the background source Cas A. Channel D has inadequate signal to noise 
to constrain the power spectrum due to very large optical depth. Model fit 
to the data of channel C (Vlsr = — 33 km s _1 ) is shown in Figure[5] 

which has been interpreted as arising from magnetohydrodynamic 
turbulence in this shell type supernova remnant. The angular scale 
of the break in the power spec trum corresponds to the shell thick- 
ness of Cas A jRov et alj|2009t) . 



4 THE OPACITY FLUCTUATION POWER SPECTRUM 

The sky brightness S(l, to) towards the extended source Cas A may 
be written as 

S(l,m) = J c (Z,m) e - T(i ' m) (4) 

where I c (l, m) is the continuum radiation from Cas A and t(1, to) 
is the optical depth for absorption by intervening H I. We expect 
I c (l, to) and t(7, m) to be independent, whereby we can write 

= (I c (l,m)I c (l',m')) x (e-^."0 e - T < , '' m ' ) ). (5) 

Both I c (l,m) and r(7,m) can be decomposed into the sum of a 
constant average value and the fluctuation around that value (as a 
function of I, to) 

I c (l,m) = I° + SI c (l,m) 

r{l,m) = T°+Sr(l,m) (6) 
and hence equation $5$ can be re-written as 

^ = [I c +Zlc) X e ( e e > ( ? ) 

where £/ c is the autocorrelation function of continuum intensity 
fluctuation. Further, we assume that the opacity fluctuation 5r is a 
Gaussian random field with zero mean, variance a 2 and autocorre- 
lation function £ T , whereby 

( e -6T(l,m) e -ST(l',m')} = exp((J 2 + ^ _ (8) 

The opacity power spectrum P T — T^t), which is our main inter- 
est in this paper, is the Fourier transform of £ T . The power spectrum 
of the brightness fluctuations for a frequency channel with H I ab- 
sorption can be written as 

Ps = H^} = e- 2T ° + ^(I° c 2 + P Ic ) ®T{e^} (9) 
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where Pj c is the power spectrum of brightness fluctuations of the 
continuum emission, i.e. the Fourier Transform of £j c . Equation ([9} 
relates Ps to the convolution of (7° 2 + Pi c ) and _F{e^ T } which 
refer to the continuum and the H I opacity respectively. It should 
be noted that while the terms Ps, Pi a and T{e^ T } are all functions 
of U, the terms 7° and e~ 2r +<Tt do not have any U dependence. 

_o I 2 

In our analysis we have treated B = e T as a free param- 

eter that determines the amplitude of Ps- The terms 7° 2 and Pi c , 
that refer to the continuum radiation, are known from Cas A's total 
continuum flux density and analysis of the brightness fluctuations 
in the line free channels (i.e. eq. [3} respectively. The H I opacity 
power spectrum P T , which is the quantity that we wish to deter- 
mine in this analysis, enters eq. ((9) through the term T{e^ T }. We 
have assumed that the opacity power spectrum is a power law 



Pr =-^{£t} = AIT. 



(10) 



Under this assumption the brightness power spectrum Pg(U) is 
completely determined by three parameters A, B and a. For each 
velocity channel that shows H I absorption we have used standard 
chi-square minimization to determine the value of these three pa- 
rameters for which our model prediction (eqs.l9landll0t best fits the 
measured Ps- The values of A and B are not of particular interest 
in the present analysis. Hence we derive only the best fit value of a 
and marginalize over A and B to estimate errors for a. 



5 RESULTS AND DISCUSSION 

Figure[3]shows the measured Ps and the best fit model of Ps for a 
single channel with H I absorption. We find that our model, which 
assumes a power law H I opacity power spectrum (eq.UOK provides 
a good fit to the measured sky brightness power spectrum Ps. The 
reduced chi-square values for the best fit model are, however, not 
close to 1.0 but are about 0.3. This is most probably because of 
an overestimation of the errors in the measured Ps. The estimated 
errors, for example, are found to be significantly more than the vari- 
ation of the continuum power spectra across line free channels. It 
implies that the very conservative and approximate error estimates 
adopted here may result in a smaller value of the reduced chi-square 
even when the model represents the observed power spectrum quite 
well. Hence, we have scaled down the errors of the observed Ps so 
that the reduced chi-square is close to 1.0. The la errors of the 
power law index correspond to a change of the reduced chi-square 
by +1 from its minimum value. 

The velocity channels with very strong absorption do not have 
adequate signal to noise to constrain the power spectrum. On the 
other hand, for a few channels with dominant continuum signal and 
very weak H I absorption, the power spectrum does not deviate sig- 
nificantly from the "continuum" power spectrum, making it hard to 
constrain the opacity spectrum. The optimal channels for this anal- 
ysis are hence those with moderate H I absorption, where there is 
still sufficient signal to noise ratio, and the total power spectrum has 
been significantly modified by the H I absorption. We note that any 
image-based method of estimating the power spectrum will also 
be constrained to the same optimal channels by the signal to noise 
considerations mentioned here. 

The best fit value of a is found to be —2.86 ± 0.30 (3cr er- 
ror) for the angular scale of 5 — 150" probed in this observation. 
Note that this error bar is based on the quadrature sum of the es- 
timated noise due to "cosmic variance" and the signal to noise ra- 
tio on each visibility. As discussed above, this is likely to be an 
overestimate of the true error. Consistent with this, the channel 
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a = -2.8 + 0.3 
V LSR = -33 km 5' 
AV = 0.4 km s" 1 




Figure 3. Intensity fluctuation power spectrum for a "thin" velocity channel 
with H I absorption (channel C marked in FigurefJJand FigureffJ- The best 
fit model power spectrum is plotted with the data points. 



to channel variation of the best fit value of a is much less than 
the estimated error of 0.30. Computing the variance over indepen- 
dent estimates of a derived for different velocity channels leads to 
a ~ —2.86 ± 0.10 (3<7 error). We find no significant difference in 
a between the velocity channels with absorption produced by gas 
from the Perseus arm and from the Local arm. The angular scales 
that our power spectral estimator is valid for corresponds to a linear 
scale of 0.07 — 2.3 pc for the Perseus arm (for an assumed distance 
of 3 kpc to the ar m). The Local arm is ~ 100 pc wide in the so- 
lar neighbourhood jPorcel etalj|l998h . So, for the Local arm, the 
upper limit of linear scale range (adopting a distance of 100 pc to 
the absorbing gas) is 0.002 - 0.07 pc (48 - 15730 au). We note 
that, for the Local arm towards Cygnus A. lDeshpande et all 12000) 
have reported a = —2.5, an index somewhat shallower than what 
we find here, or what those authors find for the Pe rseus arm and the 
Outer arm. For the Perseus arm towards Cas A, Deshpande et al. 
(2000) report a = —2.75 ± 0.25 (3<r error) for a similar range of 
linear scales. This factor of ~ 2.5 smaller uncertainty in the current 
measurement is partly due to the better sensitivity of the present ob- 
servations and partly due to using a more optimal estimator of the 
power spectrum. The value of a that we find here is also consis- 
tent with the value (~ —3) reported from the Galactic H I emission 
observations (Green 1993) probing the scales of 50 — 200 pc. 

For power law opacity fluctuations power spectrum (eq. ITOb 
with —2 a —4, the o pacity fluctuation smoothed over 
a length-scale x is predicted (Desh pandell200(]|) to have a rms 

a — 2 

a T (x) oc x~5~ which is also a power law. Using our best fit 
model we have a T (x) = 2.0 (x/A pc) ' 43 . Extrapolating this 
power law, the rms opacity variation is ~ 0.05 on scales of 100 
au which is broadly consistent with recent VLA and VLBA studies 
dFaison etalJ|l99St iDeshpande et alj|2000l ; iFaison & Go"ssl 1200 ll; 
iBrogan et al . 2005; L azio et al.l 12009). The present study of H I 
absorption, combined with this result derived from H I emission 
observations at larger linear scales indicates that there is no signifi- 
cant change in the slope of the power spectrum over 4 — 5 orders of 
magnitude in linear scale. This spectrum is significantly shallower 
than the Kolmogorov spectrum (with power law in dex of 11/3) 
expec ted from an incompressible turbulent medium ( Kolmogorov 
1 19411) . This may be because of the fact that the turbulen ce in H I is 
comp ressible and magnetohydrodynamic in nature (see lRov et al.l 
2009, for discussion on this issue). 
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If the density fluctuation is small compared to the mean den- 
sity, then for gas in pressure equilibrium, the slope of the opacity 
fluctuation power spectrum will be nearly the same as tha t of the 
density fluctuation power spectrum dDeshpande et al.l l2000). How- 
ever, fluctuations in both the density and velocity fields contribute 
to the observed H I opacity fluctuations. It can be shown that for 
the power spectrum estimated from "thick slices", i.e. those with 
velocity width larger than the turbulent velocity dispersion, the con- 
tribution is only from the de nsity fluctuations, and all ve locity in- 
formation get averaged out dLazarian &Pogosvarjl2000h . So, for 
"thick slices", the intensity fluctuation is dominantly due to density 
fluctuation and the observed power law index n is expected to be 
same as the index of the density fluctuation power spectrum. On the 
other hand, for "thin slices", the power law index of the observed 
power spectrum is n + 13/ 2, where /3 is the power law index of 
the velocity structure function. Hence, it is possible to decouple the 
density and velocity fluctuations from the power spectra derived 
from "thin" and "thick" velocity channels. lLazarian"& Pogosvan 
(2000) reported a near-Kolmogorov power law index for both den- 
sity and velocity fluctuations power spectra for the Milky Way and 
the Small Ma gellanic Cloud by applying this techni que to the H I 
emission data (Green 1993; Stanimirovic et aljl999h . For our data, 
we find that though there is a very weak trend of smaller a for 
larger channel width, the observed change of a by about —0.1 for 
a change of velocity width from ~ 0.4 to ~ 12.8 km s _1 is within 
the 1 a error bar. Since a remains unchanged for "thin slices" with 
velocity width as small as ~ 0.4 km s , i.e. much smaller than the 
typi cal value of 4.0 km s" 1 for turbulent dispersion in Galactic cold 
H I iRadhakrishn an et aljl972t) . we can constrain /3 to be 0.2 ± 0.6 
(3cr error). This is consistent with the value of /3 = 2/3 predict ed 
for turbulence in an incompressible medium I Kolmog orovll 1 94 il) . 



6 CONCLUSIONS 

In this work, we have studied the H I 21 cm opacity fluctuation 
towards Cas A. A simple but robust method of directly estimat- 
ing the power spectrum from the observed visibilities is outlined. 
In this analysis we avoid the complications of imaging the bright, 
extended continuum source, of subtracting the continuum and of 
making optical depth image for channels with H I absorption. 

We have found that the H I opacity fluctuation power spectrum 
can be modelled as a power law with an index of —2.86 ± 0.10 
(3cr error). This is consistent with earlier observational results. We 
have not found any significant difference of the power law index 
between velocity channels with absorption produced by the gas 
from the Perseus arm and the Local arm. We have also checked, 
by smoothing the visibility data for adjacent channels, if there is 
variation of the power law index with the velocity width of chan- 
nels. It is found that, within the estimation errors, the power law 
index remains constant for a wide range of velocity widths. We can 
not, however, rule out contribution of velocity fluctuations to the 
observed opacity fluctuations if the power law index of the velocity 
structure function is 0.2 ± 0.6 (3(7 error). We plan to apply, in fu- 
ture, this visibility based method to study the opacity fluctuations 
for other lines of sight in different regions of the Galaxy. 
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